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ABSTRACT. Solutions to conservation laws satisfy the monotonicity property: 
the number of local extrema is a non-increasing function of time, and local maxi- 
mum/minimum values decrease/increase monotonically in time. This paper inves- 
tigates this property from a numerical standpoint. We introduce a class of fully 
discrete in space and time, high order accurate, difference schemes, called gener- 
alized monotone schemes. Convergence toward the entropy solution is proven via 
a new technique of proof, assuming that the initial data has a finite number of 
extremum values only, and the flux-function is strictly convex. We define discrete 
paths of extrema by tracking local extremum values in the approximate solution. In 
the course of the analysis we establish the pointwise convergence of the trace of the 
solution along a path of extremum. As a corollary, we obtain a proof of convergence 
for a MUSCL-type scheme being second order accurate away from sonic points and 
extrema. 
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1. Introduction 

This paper deals with entropy solutions to the Cauchy problem for a one- 
dimensional scalar conservation law: 

d t u + d x f(u) = 0, u(x,t) e R, t > 0, x e R, (1.1) 

u(x, 0) = «o(a;), x E R, (1-2) 

where the flux / : R — > R is a given function of class C 2 and the initial data uq 
belongs to the space BV(R) of all integrable functions of bounded total variation. 
For the main result of this paper, we assume that 

/ is a strictly convex function (1.3) 

and 

uq has a locally finite number of extrema. (1-4) 

Solutions to conservation laws are generally discontinuous and an entropy criterion 
is necessary to single out a unique solution. We refer the reader to Lax [20, 21] for 
background on nonlinear hyperbolic equations and the entropy criterion. 

As is well-known [19, 37], problem (1.1)-(1.2) admits a unique entropy solution 
u in the space L°° (R + , BV(R)) and Lip(i? + ; L 1 (.#)). (This result holds without 
the restriction (1.3)-(1.4) and for multidimensional equations as well.) To select 
the solution one can use the distributional entropy inequality 

d t U(u) + d x F(u)<0, (1.5) 

where the (Lipschitz continuous) function (£/", F) : R — > R 2 is a convex entropy- 
entropy flux pair, i.e. U is a convex function and F' = U'f. One can also use the 
Lax shock admissibility inequality 

u(x-,t) > u(x+,t) (1.6) 

for all x and t. Since / is convex and u has bounded variation, a single entropy 
U in (1.5) is sufficient to ensure uniqueness, and (1.6) is equivalent to (1.5). At 
the discrete level, however, conditions (1.5) and (1.6) leads two drastically different 
notions of consistency with the entropy criterion for a difference scheme. In the 
present paper we will be using both conditions. Indeed in some regions of the (x,t) 
plane it is easier to use a discrete version of (1.5), while in other regions (1.6) is 
more adapted. 

We are interested in conservative discretizations of problem (1.1)-(1.2) in the 
sense of Lax and Wendroff [22]. The monotone schemes and, more generally, the 
E-schemes are large classes of schemes (including the Godunov and Lax-Friedrichs 
schemes) for which a convergence analysis is available. The main point is that 
monotone schemes satisfy a discrete version of the entropy inequality (1.5) (see 
(2.9) below). However they turn out to be first order accurate only, and so have 
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the disadvantage to introduce a large amount of numerical viscosity that spreads 
the discontinuities over a large number of computational cells. 

The proof of convergence of the monotone schemes and E-schemes is based on 
Helly's and Lax-Wendroff 's theorems. See [3, 6, 8, 11, 16, 26, 31, 35, 23] and the 
references therein. The result holds even for multidimensional equations and/or 
when irregular (non-Cartesian) meshes are used. 

To get high-order accurate approximations, it is natural to proceed from ana- 
lytical properties satisfied by the entropy solutions to (1.1), formulate them at the 
discrete level, and so design large classes of high-order difference schemes. One 
central contribution in this direction is due to Harten [13, 14], who introduced the 
concept of TVD schemes, for "Total Variation Diminishing". Harten shows that 
conservative, consistent, TVD schemes necessarily converge to a weak solution to 
(1.1). Moreover such schemes possess sharp numerical shock profiles with no spu- 
rious oscillations. However, Harten's notion of TVD scheme is weaker than the 
notion of monotone scheme, and a TVD need not converge to the entropy solution. 
The aim of this paper is precisely to single out a subclass of TVD schemes, refining 
Harten's notion, that are both high-order order accurate and entropy satisfying, 
cf. Definition 2.2 below. 

A very large literature is available on the actual design of second-order shock- 
capturing schemes. One approach to upgrading a first order scheme was proposed 
by van Leer [23, 24]: the MUSCL scheme extends the Godunov scheme by replacing 
the piecewise constant approximation in the latter with a piecewise affine approxi- 
mation. The heart of the matter is to avoid the formation of spurious oscillations 
near discontinuities. This is achieved by van Leer via the so-called min-mod limitor. 

Other classes of high order schemes have been built from the maximum principle 
and a monotony condition: see the classes of ENO and UNO schemes proposed 
by Harten, Engquist, Osher, and Chakravarthy [15]. In this paper we concentrate 
attention on the MUSCL scheme, but our main convergence theorem, Theorem 2.3 
below, should also be applicable to other schemes. 

After Osher's pioneering result [32] and the systematic study of high-order schemes 
by Osher and Tadmor [33] (which however required a technical condition on the 
slopes of the affine reconstructions), Lions and Sougadinis [27, 28] and, indepen- 
dently, Yang [38] established the convergence of the MUSCL scheme. Both proofs 
apply to an arbitrary BV initial data (so (1.4) is not assumed); the flux-function is 
assumed to be convex in [38] and strictly convex in [28] . In both papers significantly 
new techniques of proof are introduced by the authors. In [38], Yang develops a 
method for tracking local extremum values. In [28] , Lions and Sougadinis elegantly 
re-formulate the MUSCL scheme at the level of the Hamilton-Jacobi equation asso- 
ciated with (1.1) and rely on Crandall-Lions' theory of viscosity solutions for such 
equations. A large class of difference approximations is treated in [28]. Techniques 
in both papers are restricted in an essential way to semi-discrete schemes, in which 
the time variable is kept continuous. In [25] the present authors announced a proof 
of the convergence of a class of fully discrete schemes that include the MUSCL 
scheme. Independently, Yang [39] also extended his approach to a large class of 
fully discrete methods. 
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Several other works deal with the convergence of the van Leer's scheme or vari- 
ants of it. A discrete version of inequality (1.6) was established by Brenier and 
Osher [2] and Goodman and LeVeque [12], the latter dealing with both first and 
second order methods. Nessyahu, Tadmor, and Tamir [30] establish both the con- 
vergence and error estimates for a variety of Godunov-type schemes. Various Ap- 
proaches to upgrading Lax-Friedrichs scheme are actively developed by Tadmor and 
co-authors; see, for instance, Nessyahu and Tadmor [29]. An extensive discussion of 
the discretization of the entropy inequality (1.5) is found in Osher and Tadmor [33] 
and the many references cited theirein. See also Bouchut, Bourdarias and Perthame 
[1], and Coquel and LeFloch [4]. 

The objective of the present paper is to provide a framework to prove the conver- 
gence of high order accurate and fully discrete difference schemes. We will strongly 
rely on a property shared by all entropy solutions to (1.1): the monotonicity prop- 
erty. Given an arbitrary entropy solution u, the number of extrema in u(t) is a 
non-increasing function of t, and local maximum/minimum values decrease/increase 
monotonically in time. (See the Appendix for rigorous statements.) This property 
was studied first by Harten [13] from the numerical standpoint, in order to arrive 
at his notion of TVD scheme. It was also essential in [28, 38] and in Tadmor [36]. 

Observe that monotone or TVD schemes do not necessarily satisfy the mono- 
tonicity property. For instance the Lax-Friedrichs scheme may increase the num- 
ber of extremum values! This motivates the introduction of a subclass of TVD 
schemes, guaranteeing this property together with the high order of accuracy. The 
scheme then closely mimics an important property of the solutions to the con- 
tinuous equation. Further requirement is necessary to ensure that the scheme is 
entropy-satisfying. 

Based on the monotonicity property, we thus introduce the notion of generalized 
monotone scheme (Definition 2.2) characterized by three conditions: maximum 
principle, entropy consistency, and local behavior at extrema. The (first order) 
Godunov method and the (second order) MUSCL method are prototype examples. 

Our first requirement is a (strong) version of the local maximum principle, which 
will allow us to show a discrete analogue of the monotonicity property. In particular 
our condition prevents the formation of spurious numerical oscillations and the 
scheme is TVD in the sense of Harten. 

The second requirement is motivated by the following observation made by Os- 
her for semi-discrete schemes [31]: in the non-decreasing parts of the approximate 
solutions, it is possible to simultaneously achieve second order accuracy and the ex- 
istence of one cell entropy inequality -i.e. a discrete analogue of (1.5). Our second 
condition therefore requires a cell entropy inequality in the non- decreasing regions. 

Finally, in order to prevent cusp-like behavior near extrema, which may lead 
to entropy violating discontinuities, the scheme should be well-behaved near local 
extrema. To this end, we introduce a condition referred to as the quadratic decay 
property at local extrema. It reflects the nonlinear behavior of the numerical flux- 
function near extrema, and assumption (1.3) again is essential. For simplicity, it 
is also assumed that the scheme reduces to a three points, first order scheme at 
extrema. 
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The rest of this section is devoted to comments upon the proof of our main result 
that any generalized monotone scheme converges to the entropy solution of (1.1)- 
(1.2) provided (1.3)-(1.4) holds. Our approach was driven by Yang's paper [38] on 
semi-discrete schemes. However, our technical arguments differ substantially from 
the ones in [38]. In particular we restrict attention to initial data having a locally 
finite number of extrema, making the tracking of paths of extremum values almost 
a trivial matter. The main part of our proof is studying the convergence of the 
traces of the approximate solution along it. We make use of the quadratic decay 
property above to exclude the formation of "cusp" near extremum values, which 
could lead to entropy- violating shock. This is the main contribution of this present 
paper. We do not believe that the extension of our proof to arbitrary BV data is 
straightforward. The result we obtain seems satisfactory however since condition 

(1.4) covers "generic" initial data. 

The proof distinguishes between the non-increasing parts and non-decreasing 
parts of the solution, and is based on several observations as follows. We use the 
notation u h for the approximate solutions, v for the limiting solution, and h for the 
mesh size. 

First of all, the strong maximum principle ensures that the u h 's are total vari- 
ation diminishing in time, so of uniformly bounded total variation. The strong 
maximum principle enables to easily define a discrete path and each time step the 
discrete path move at most one grid points. As a consequence, the discrete path 
is Lipschitz continuous. (This is a major difference between the present paper and 
Yang's paper in which the construction of the path and the limiting paths are much 
more involved.) 

By Helly's Theorem, the scheme converges in the strong L 1 topology to a limit- 
ing function, say v, which according to Lax-Wendroff 's theorem is a weak solution 
to (1.1)-(1.2). It remains to prove that v is the entropy solution. We show that 
the strong maximum principle in fact implies a discrete analogue of the monotonic- 
ity property. Relying on (1.4), we construct a (locally) finite family of Lipschitz 
continuous paths in the plane by tracking the local extrema in u h . The paths are 
shown to converge in the uniform topology to limiting curves. 

Next we make the following two observations. On one hand, in a non-increasing 
region for u h , the function v is also non-increasing, and so can only admit non- 
increasing jumps. Thus v satisfies the Lax shock admissibility inequality (1.6) in 
the non-increasing regions. On the other hand, a discrete cell entropy inequality, 
by assumption, holds in the non-decreasing regions of u h . So v satisfies (1.5) in the 
non-decreasing regions. 

It remains to prove that v has only non-increasing jumps along any path of 
extrema. This is the most interesting part of the proof. Let ifj h be an approximate 
path of extrema, and let ip be its uniform limit. The path i/j is the boundary 
separating two regions where the analysis in the paragraph above applies. Indeed 

(1.5) holds in the side where v is non-increasing and (1.6) holds in the side where 
v is non-decreasing. A specific proof must be provided to determine the behavior 
of v along the path. We analyze the entropy production in a small region of the 
plane limited on one side by the path ip h . In the course of this proof we derive a 
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uniform bound for the time integral of the local oscillation in space along the path 
ip h , which is a direct consequence of the quadratic decay property mentioned above. 
For the sake of simplicity, we use here the assumption that the scheme reduces to 
a three points, first order scheme at extrema. 

Note that assumption (1.3) is not used in the construction of the extremum 
paths, but is essential in the convergence analysis which strongly relies upon (1.6), 
only valid for convex fluxes. 

Our analysis shows that, for a class of difference schemes, certain approximate 
generalized characteristics -those issued from an extremum point of the initial data- 
can be constructed for scalar conservation laws with convex flux. Constructing 
approximate generalized characteristics issued from an arbitrary point remains a 
challenging open problem. Recall that, for the random choice scheme, Glimm- 
Lax [10] did obtain a general theory of approximate generalized characteristics 
(applicable to systems, as well). 

The organization of this paper is as follows. In Section 2, we define the class of 
generalized monotone schemes and state the main result of convergence, cf. Theo- 
rem 2.3. Section 3 contains the proof of the main result. In Section 4, we apply 
Theorem 2.3 to the MUSCL scheme and provide some additional remarks. 

2. Generalized Monotone Schemes 

This section introduces a class of TVD schemes which are built to closely mimic 
an essential property of the entropy solutions to (1.1), i.e. the monotonicity prop- 
erty: the number of local extrema is a non-increasing function of time, and local 
maximum/minimum values decrease/increase monotonically in time. See Appen- 
dix for a precise statement. We investigate here this property from a numerical 
standpoint. Monotone and TVD schemes actually do not necessarily satisfy this 
property (see Section 4 for an example), and a more restricted class, the generalized 
monotone schemes, is natural. 

We consider a (2k + l)-point difference scheme in conservation form for the 
approximation of (1.1)-(1.2): 

u] +1 = u]-\(g? +1/2 -g?_ 1/2 ), n>0,jeZ, (2.1) 

where we use the notation g^ + i/ 2 = 9{ uT j-k+i^ • • • > uT j+k) an< ^ A = r//i for the ratio 
of the time-increment r by the space-increment h. We set t n = nr, Xj = jh, and 
Xj + i/ 2 = (j + 1/2) h. The value u™ presumably is an approximation of the exact 
solution at the point (xj, t n ). As is usual, the numerical flux g : R 2k — > R is assumed 
to be locally Lipschitz continuous and consistent with /, i.e. g(v, . . . ,v) = f(v) for 
all v. Note that g may depend on A. For definiteness, we set 

u °j = T u o(v) <*y. (2-2) 

This is sufficient for second order accuracy. For higher orders, one should use a 
Runge-Kutta time-step method (Shu [34]). We also define the piecewise constant 
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function u h : R x R + — > i? by 

u h (x,t) = it?, t n <t < t n+1 , Xj_i/2 < x < x j+1 / 2 . (2.3) 

By construction -u' 1 is a right continuous function. For simplicity we assume the 
following CFL restriction: 

Asup|/»|<l/4. (2.4) 

V 

Several of the properties below would still hold if, in (2.4), one replaces 1/4 by 1 
although the proofs then become less clear geometrically. 

The (first order) Godunov scheme is based on exact solutions to (1.1), and is a 
good prototypical example to lead us to defining a class of (high order) schemes 
consistent with the monotonicity property. 

A main ingredient is the Riemann problem. Given two states v and w, we define 
R(.;v,w) to be the entropy solution to (1.1)-(1.2) with, here, 

uq(x) = v if x < 0, w if x > 0. 

As is well-known, R(-;v,w) depends on the self similarity variable x/t only, and is 
given by a closed formula. If v < w, R is a rarefaction wave and, if v > w, a shock 
wave. More important 

x 

R( — ;v, w) is a monotone function connecting v to w. (2-5) 

t 

In Godunov scheme, one solves Riemann problems and, at each time level, one 
projects the solution on the space of piecewise constant functions. If u h (t n ) is 
known, let u(x,t) for t > t n be the entropy solution to (1.1) assuming the Cauchy 
data 

u(t n +) = u h (t n +). 

Since u (t n +) is a piecewise constant function, u is obtained explicitly by glue- 
ing together the Riemann solutions R(.; it™, w" +1 ). In view of (2.4), there is no 
interaction between two nearby solutions, at least for t G (t n ,t n +i). Set 

u (x,t n+l +) = t u(y,t n+1 )dy, x j _ 1/2 <x< x j+1/2 . 

Using the conservative form of (1.1), the scheme can be written in the form (2.1) 
with k = 1 and g = go given by 

go(v,w) = f(R(0+;v,w)) for all v and w. (2.6) 

The proofs of (2.7)-(2.10) stated below are classical matter; e.g. [3, 6, 16, 26]. 
On one hand, one can think of the Godunov scheme geometrically as a two-step 
method: a marching step based on exact (Riemann) solutions and an L 2 projection 
step. In view of (2.5) and the monotonicity property of the L 2 projection, one 
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easily have a simple geometrical proof of the properties listed below. On the other 
hand, an algebraic approach is based on the explicit formula deduced from (2.6). 

The Godunov scheme is monotone: the function go is non-decreasing with re- 
spect to its first argument, and non- increasing with respect to its second one. This 
property implies that the scheme is monotonicity preserving, i.e., 

if w™, Uj 1+1 , . . . , Uj 2 is a non-increasing (resp. non-decreasing) sequence 
for some indices ji < j 2 , so is the sequence u^+i, w ? 1 "+2) • • • > u j^~-i- 

The Godunov scheme satisfies the local maximum principle, i.e., 

min u], u] +1 ) < u] +1 < max (u^, u], u] +1 ) (2.8) 

for all n > and j G Z. In fact (2.7) and (2.8) are shared by both steps in the 
Godunov scheme. It will convenient to us to rewrite (2.8) in term of the jumps of 
u h at the endpoints of a cell: 

min (u] +1 - v$, u n 3 _ x - u], 0) < u] +1 - u] < max (u] +1 - u?, u]_ x - v$, 0) . (2.8') 

Any monotone scheme -in particular the Godunov scheme- satisfies a discrete 
analogue of (1.5) for every convex entropy pair (U,F): 

U(u] +1 )-U(u])-X{G] +1/2 -G]_ 1/2 )<0, n>0,j€Z. (2.9) 

In (2.9), G» +1/2 = G(«»_ 

.fc+i) • • • ? and G is a numerical entropy flux consis- 

tent with F, that is G(v, v . . . , v) = F(v) for all v. 

Finally concerning the local behavior of u h in the neighborhood of local extrema, 
it is known that, say for local maximum, 

u n+l < w n_ ( 2 

A similar property holds for local minima. 

In fact the classical properties (2.8) and (2.10) can be improved as follows. 

Proposition 2.1. Under assumptions (1.3) and (2.4), the Godunov scheme satis- 
fies the following two properties: 

(1) the strong local maximum principle; 

l - min - u» u]_, - U y , 0) < - < 1 max - t£, tx^ - u» 0) 

(2.11) 

(2) and £/ie quadratic decay property at local extrema, that is e.g. for a maxima: 
If it™ is a local maximum value , < — amin((«" — it™ ±1 ) 2 ). (2.12) 
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with a = Ainf/72. 



□ 



The proof of (2.11) is straightforward from a geometrical standpoint. If also 
follows from Proposition 4.1 established later in Section 4. Note that the coefficient 
1/4 in (2.4) is essential for (2.11) to hold. Observe that (2.11) is stronger than 
(2.8)-(2.8') and controls the time-increment (i.e., u™ +1 — -u™') of the solution in the 
cell j in term of the jumps at the endpoints: the values evolves "slowly" as t n 
increases. As we shall see, this property implies that the scheme satisfies a discrete 
analogue of the monotonicity property. 

Estimate (2.12) is stronger than (2.10) and shows that the decrease/increase of 
a maximum/minimum is controlled by the quadratic oscillation of u nearby this 
extremum. It is a truly nonlinear property of the Godunov flux. It will be used 
below to prove that cusp can not form near extremum points. For convenience, the 
proof of (2.12) is postponed to Section 4, where second-order approximations are 
treated as well. 

In [12], Goodman and LeVeque derive for the Godunov method a discrete version 
of the Oleinik entropy inequality. In particular, this shows that the Godunov scheme 
spreads rarefaction waves at the correct rate. Our estimate (2.12) is, at least in 
spririt, similar to this spreading estimate, and expresses the spreading of extremum 
values. 

We are now ready to introduce a class of high-order schemes based on the prop- 
erties derived in Proposition 2.1. 

Definition 2.2. The scheme (2.1) is said to be a generalized monotone scheme if 
any sequence [u™ } generated by (2.1) satisfies the following three conditions: 

(1) the strong local maximum principle (2.11), 

(2) the cell entropy inequality (2.9) for one strictly convex pair (U, F) in any 
non- decreasing region, including local extrema, 

(3) the quadratic decay property at local extrema (2.12) for some constant a > 0. 

It is also assumed that the numerical flux and the numerical entropy flux are es- 
sentially two-point functions at local extrema. □ 

According Proposition 2.1, the (first order) Godunov scheme belongs to the class 
described in Definition 2.2. Section 4 will show that there exist high order accurate 
schemes satisfying the conditions in Definition 2.2. Our main convergence result is: 

Theorem 2.3. Let (2.1) be a generalized monotone scheme. Assume that assump- 
tions (1.3)-(1.4) hold together with (2.4). Then the scheme (2.1) 

(1) is L°° stable, i.e., 



inf uf < 



n+1 



< supw[\ 
lez 



n > 0, j e Z, 



(2.13) 



(2) is total variation diminishing, i.e., 




n > 0, 



(2.14) 



jez 
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(3) and converges in the Lf oc strong topology for all p G [1, oo) to the entropy 
solution of (1.1) -(1.2). □ 

The proof of Theorem 2.3 is given in Section 3. 

Theorem 2.3 is satisfactory for a practical standpoint. Suppose that uq is an 
arbitrary BV function, and we wish to compute an approximation to the solution u 
of (1.1)-(1.2) of order e > in the L 1 norm. Let us determine first an approximation 
of «o, say wo, e , that has a finite number of local extrema and such that 

IK, £ - wo II l 1 (J?) < £■ 

Applying a generalized monotone scheme to the initial condition uq jS yields an 
approximate solution that, in view of Theorem 2.3, satisfies 

||«e - u e \\ L i {R) < o(h) < e 

for h is small enough, where u e is the entropy solution associated with the initial 
condition w , e . Since the semi-group of solutions associated with (1.1) satisfies the 
L 1 contraction property, one has 

\\u £ ~ u\\ L i( R ) < |K,e - U \\ L i( R ) < 

and therefore 

\\ u e ~ u \\l^(R) < - U s \\ L i( R ) + \\u e ~ u\\ L i(R) 

<2e. 

3. Convergence Analysis 

The proof of Theorem 2.3 is decomposed into several lemmas, Lemmas 3.1-3.12. 
For the whole of this section, we assume that the hypotheses made in Theorem 2.3 
are satisfied. 

We introduce first some notation and terminology. We call u™ a local maximum 
or a local minimum if there exist two indices an d j* with j* < j < j* such that 

u£ = = --- = u^> max(^_ 1 ,w^ +1 ) 

or 

ul = ul +1 = ... = «»< min^.!,^)- 

In such a case, there is no need to distinguish between the extrema , 
. . . , u™*. Based on the strong maximum principle (2.11), we show in Lemmas 3.1 
and 3.2 that the scheme satisfies a discrete form of the monotonicity property. We 
construct a family of paths in the (x, t)-plane by tracing in time the points where 
the approximate solution u (t) achieves its local extremum values. One difficulty 
is proving that the interaction of two (or more) paths does not create new paths, 
so the total number of paths at any given time remains less or equal to the initial 
number of local extrema in uq. In passing we observe that an extremum point 
moves one grid point at each time-step, at most. 
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Lemma 3.1. For some j* < j* , suppose that the sequences _ 3 , Uj t _ 2) Uj.-n 
and u™*, u™* +1 , u™* +2 > w™* +3 are two monotone sequences, no specific assump- 
tion being made on the values u™ , j* < j < j* . Then the number v' of extrema in 
the sequence 

S n+1 := («? + ) i ,_ 2 < J -<j.+2 
is less or equal to the number v of extrema in 

S n ■= (uj ) j t -2<j<j*+2• 
When v' > 1, there exists a one-to-one correspondence between v' local extrema 
of S n and the v' local extrema of S n+ i with the following property: if a maxi- 
mum/minimum Uj is associated with a maximum/minimum u™, +1 , then 

< 1 and < u], resp. > u]. (3.1) 

□ 

Proof. We distinguish between various cases depending on the number of local 
extrema in the sequence S n and construct the one-to-one correspondence. 

If S n has no local extremum, for instance is non-decreasing, then S n +i is also 
non-decreasing. This indeed follows from inequalities (2.11) which reduce in this 
case to 



Consider next the case that S n has exactly one local extremum, say a local 
maximum at some uf. The same argument as above shows that the sequences 
. ^ . , and {u" --1-1 ), , , . are non-decreasing and non-increasing re- 

spectively. Therefore we only need to exclude the case that both > up +1 

and itj^ 1 > u™ + , which would violate the monotonicity property since Sn-\-i in 
this case would have two local maximum and one local minimum, so two new ex- 
trema. Indeed assume that the latter would hold, then using (2.11) at the points 
1 — 1,1, and I + 1 gives us 



u 



n+l 



1 



«P + - min («P ±1 - uP) <uP +1 , 



and 



u 



n+l 



u? +1 



which are incompatible with the inequalities upj^ 1 > u™ +1 and u^ 1 > u, z 

Consider now the case that 5 n has two local extrema, say one local maximum 
at / and one local minimum at m with / < m. We distinguish between three cases: 
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If I < m — 2, then the two extrema can not "interact" and the arguments before 
show that the solution at time £ n +i has the same properties. 

If / = m — 2, the two extrema can interact. Using (2.11) at each point j = 
/ — !,••• , / + 3, one gets 



u 



n+l 
l- 



1 



1 <o K-l+O' 



1 



u? + £ min («? ±1 - «?) <<+\ 



and 



and 



it 



n+l 



1 



i 1 <o(«r+i+«r). 



i 

2 

1 



< + + 1 1 <^(< + l+< +2 ): 



U 



and 



r +2 + g min i u ?+2±i - <+ 2 ) < 



2 

It is not hard to see that these inequalities imply that S n+ i 

(1) either has one maximum at either j = 1 — 1,1,1 + 1 and a minimum at 
j = 1 + 1,1 + 2,1 + 3, 

(2) or is non-decreasing. 

In the first case, we achieve the property we wanted. In the second case, there is 
no extremum at the time t n +i- 

This analysis can be extended to the case that several extrema can "interact" ; we 
omit the details. Property (3.1) is a consequence of condition (2.4): an extremum 
point can only move up to one grid point at each time step. □ 

Consider the initial condition uq and its approximation u h (0) defined by L 2 
projection, cf. (2.2). Locate the minimum and maximum values in the initial data 
uq. For h much smaller than the minimal distance between two extrema, ^(O) has 
the same number of extrema as uq and the same increasing/decreasing behavior as 
uo. Indeed, there exist indices Jg(0) for q in a set of consecutive integers E{uq) 
depending on uq but not on h, such that 

is non-decreasing for J2p(0) < j < ^2 P +i(0) 
is non-increasing for J 2p _i(0) < j < ^(0)- 

Those indices are not uniquely determined in the case that uq is constant on an 
interval associated with a local extremum. Since uq has a locally finite number 
of local extrema, there exists a partition of Z into intervals (j*, j*) in which the 
hypothesis of Lemma 3.1 holds. It is an easy matter to use the one-to-one corre- 
spondence in Lemma 3.1 and trace forward in time up to time t\ = r the locations 
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of the extrema in . At each time level a (possibly new) partition of Z is 

considered and Lemma 3.1 is used again. Indeed the values Jq(n + 1) in Lemma 
3.2 below are defined from the J^(n)'s according to the one-to-one correspondence 
established in Lemma 3.1. Finally piecewise affine and continuous paths are ob- 
tained by connecting together the points of local extrema. It may happen that 
the number of extrema decreases from time t n to £ n +i- I n such a case, one path, 
at least, can no longer be further extended in time and so, for that purpose, we 
introduce a "stopping time", denoted by = t n . 
The following lemma is established. 

Lemma 3.2. There exist continuous and piecewise affine curves ip^ : [0,T^] — > R 
forq G E(uo), passing through the mesh points (xjh^ n ^t n ) and having the following 
properties: 

^qit) = x J£(n) H ~{ X J%{n+l) ~ X J%(n))i * ^ [t n ,t n+1 ], (3.3) 

for each n = 0, 1, 2, . . . , iVj with T% = N^t < oo, 

V# < V'J+n \%J%(n) - xjh (n+ i)| < /i; (3.4) 



t/iere zs on/y a /inite number (uniformly bounded w.r.t. h) of curves ip 
on each compact set 

and 



q 



(3.5) 



(3.6) 



x G (ip2p(tn)iip2p+i(tn)) l—> u h (x,t n ) is non- decreasing, 
x G (V , 2p-i(*n)>V , 2p(*n)) ^ w (x, t n ) non-increasing . 
Furthermore, the functions : [0, T^ 1 ] — > i? defined by 

w h q {t) = u n w for t n <t<t n+1 (3.7) 

are non- decreasing if q is even, and non-increasing if q is odd. □ 

Remark 3.3. 1) The definition (3.3) is not essential. All the results below still 
hold if is replaced by any (uniformly) Lipschitz continuous curve passing through 
the mesh points (xjh^, t n ). As a matter of fact, it is an open problem to show 
the strong convergence of the derivatives of approximate paths. By comparison, for 
the approximate solutions built by the random choice scheme, Glimm and Lax [10] 
prove the a.e. convergence of the first order derivatives of the paths. 
2) Introducing the stopping times T q is necessary. At those times, certain paths 
cross each other and their extension in time is not well-defined. For instance a 
path of maximum and a path of minimum can cross and "cancel out". The case 
of exact solutions (Cf. the appendix) is simpler in this respect: the paths can be 
defined to be characteristic curves for all times, even when they are no longer paths 
of extrema. 
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3) It is not interesting to trace the minimal (or maximal) paths of extrema in the 
approximate solution. Such paths would not converge to the paths obtained in the 
continuous case. □ 

By construction, cf. (3.4), a path "jumps" up to one grid point at each time-step. 
So the slope of a path remains uniformly bounded by 1/A and the curves ip^ are 

bounded in the W 7 ^, 00 norm, uniformly with respect to h and q. On the other hand, 
Lemma 3.3 implies that the scheme is TVD so TV(u h (t n )) is uniformly bounded. 
We thus conclude that the approximate solutions and the paths of extrema are 
strongly convergent, as stated in the following Lemmas 3.4 and 3.5. For simplicity, 
we keep the same notation for a sequence and a subsequence. 

Lemma 3.4. There exist times T q G [0, +oo] and Lipschitz continuous curves ip q : 
(0, T q ) -> R such that 

Ty -*T q as h -> 0, (3.8) 

and 

~^ uniformly on each compact subset of C°([0, T q )). (3.9) 

□ 

Lemma 3.5. The sequence u h satisfies estimates (2.13) -(2.14), and so is uniformly 
stable in the L°°([0, oo), BV(R)) and Lip([0, oo), ^(R)) norms. There exists a 
function v in the same spaces such that 

u h (x,t) — > v(x,t) for all times t > and almost every x G R, (3.10) 

and there exist functions w q in BV((0,T q ),R) such that 

Wq — > w q almost everywhere on (0,T 9 ) (3.11) 

for all q G E{uq). □ 

The convergence results (3.10) and (3.11) hold in particular at each point of 
continuity of v(t) and w q , respectively. Introduce now the following three sets, 
which provide us with a partition of the (x, £)-plane into increasing/decreasing 
regions for v: 

Q^V) = {(x,t) I 1p 2pi (t) < X < 1p 2p2 + 1 (t), pi <p 2 ,t < T 2pi , t < T 2p3 + i, 

and t > T q , for all q = 2pi + 1, 2p 2 }, 
Q, 2 (v) = {(x,t) /ip 2pi -i(t) < x < ip2 P2 (t), pi <p 2 ,t< T 2pi _i, t < T 2p2 , 

and t > T qi for all q = 2pi, 2p 2 — l}, 

^3(^) = Closure { (ip q (t) , t) , for all relevant values of t and q}. 

The set 03(f), by construction, contains all of the curves ip q including their end 
points. The sets fl±(v) and H 2 (f) are open and contain regions limited by curves 
in 03(f). These definitions take into account the fact that the path need not be 
defined for all times. Observe also that an arbitrary point in 3 (i> ) need not be a 
point of extremum value for v. The decomposition under consideration is not quite 
the obvious partition of the (x, t)-plane into regions of monotonicity for v. Strictly 
speaking, the sets £lj(v) may not be determined from the function v alone. 
Using Lemmas 3.2, 3.4, and 3.5, we immediately check that: 
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Lemma 3.6. The limiting functions satisfy the properties: 

u/nj („)(£) is non- decreasing in each subcomponent of Cli(v), 
v/n 2 ( v )(t) is non-increasing in each subcomponent of 2 (f), 



(3.12) 



and 



w q is non- decreasing if q is even and non-increasing if q is odd. (3.13) 

□ 

Since the scheme is consistent, conservative, and converges in the L 1 strong 
norm, we can pass to the limit in (2.1). It follows that v is a weak solution to (1.1). 
We note that, in the set Qi(v), the functions u h and, so v, are non-increasing. 
The Lax shock inequality holds for both u h and v. On the other hand, the cell 
entropy inequality (2.9) holds for u h in the non-decreasing regions, i.e., in Q,2{v). 
The passage to the limit in (2.9) is a classical matter. 

Lemma 3.7. The function v is a weak solution to equation (1.1) and satisfies 

v(x—,t) > v(x+,t) in the set Qi(v) (3-14) 

and 

d t U(v) + d x F(v) < in the set Sl 2 (v). (3.15) 

□ 

The rest of this section is devoted to proving that the Lax shock inequality holds 
along the paths ip q which we will attain in Lemma 3.10. In a first stage, we prove: 

Lemma 3.8. Along each path of extremum ifj q and for almost every t e (0,T q ), 
one of the followings hold: 

w q (t) = v(ip q (t)-,t) or w q (t) = v(ip q (t)+,t). (3.16) 

□ 

Roughly speaking (3.16) means that that no cusp-like layer can form in the 
scheme nearby local extrema. The idea of the proof of Lemma 3.8 is as follows: 
we are going to integrate the discrete form of the conservation law (2.1) on a 
(small) domain limited on one side by an approximate path of extremum, then we 
shall integrate by parts and pass to the limit as h — > 0. Finally, we shall let the 
domain shrink and reduce to the path itself. To determine the limits of the relevant 
boundary terms as h — > 0, we have to justify the passage to the limit in particular 
in the numerical fluxes evaluated along the approximate path. Lemma 3.9 below 
provides us with an a priori estimate for the oscillation of u h along the path, which 
follows from the quadratic decay property (2.12). 
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Lemma 3.9. Along a path of extremum values ip^, we have 

n=ru r 

(3 ^2 mm (l U j£(n)±l -U j£(n)l> \ U \(n)±l- U \(n)\ 2 ) ^ U (n_ ) ~ U (n+) + l ' ( 3 ' 17 ) 



n—n- 



for all < n_ < n + < , where 

/? = min(a,l/2). (3.18) 

□ 

Proof of Lemma 3.8. We will prove that, for almost every t in (0,T 9 ), the 
following three Rankine-Hugoniot like relations hold: 

' /0g -(t)(^(^ g (t)+,t)-^(^(t)-,t))+/(^(^(t)+),t)-/( V (^ g (t)-,t)) = 0, (3.19) 



dt 

(«(lMf)±. *) " «\r (*)) + /(^(*)±, *)) - /KW) = 0. (3.20) 

Since there is only one non-trivial pair of values that achieves a Rankine-Hugoniot 
relation for a scalar conservation law with a strictly convex flux and a given shock 
speed dip q (t)/dt, the desired conclusion (3.16) follows immediately from (3.19)- 
(3.20). 

Observe that (3.19) is nothing but the standard Rankine-Hugoniot relation since 
v is a weak solution to (1.1) and tf) q is Lipschitz continuous. For definiteness we 
prove (3.20) in the case of the "+" sign. The proof of (3.20) with "— " sign is 
entirely similar. (Actually only one of the two relations in (3.20) suffice for the 
present proof.) 

Let 9(x, t) be a test-function having its support included in a neighborhood of 
the curve ip q and included in the strip Rx (0, T q ). So for h small enough the support 
of 9 is included into R x (0, T^) and all the quantities to be considered below make 
sense. Let us set 9™ = 9(xj,t n ). To make use of estimate (3.17), it is necessary 

to define a "shifted" path ipq, to be used instead of ip q . So we consider the set of 
indices 

^ = {(j,n)/j>j;(n) + e h q (n)}, 

where e q (n) = (respectively e q (n) = 1) if J q (n) — 1 (resp. J q (n) + 1) achieves 
the minimum in the left hand side of (3.17). A shifted path is defined by 

^q(t) = x J*(n)+e%(n) H ^ ( X J£ (n+l)+e£ (n+1) ~ x J q h (n)+eJ(n)) ) t E [t n ,t n+1 ]. 

Introducing the shifts (n) does not modify the convergence properties of the path. 
It is not hard to see, using solely the fact that the path is uniformly bounded in 
Lipschitz norm, that as h — > 

i) h Q -> i) q W 1 ' 00 weak- ★ . (3.21) 
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We also set 

= U ^(n)+e^n) for t n < t < t n+1 . 

Using (3.17) and (3.7) of w^, it is checked that 

Wq -> w 9 L 1 strongly. (3.22) 

Consider 

I h (B) = ~ E K +1 -^ + A (^ + i/ 2 -^-i/ 2 ))^ = 0, (3.23) 

(n,j)£V% 

which vanishes identically in view of (2.1). Using summation by parts gives 

I h {9) = - £ - + E^(»He5(»)-l/2^(») T 

+ E Ui +1 W +1 -0?)b + 9i + i/2W + i-Wf)'r ( 3 - 24 ) 

= li(6) + lS(6) + lj(6). 

The passage to the limit in 7^ (0) is an easy matter, since it has the classical form 
met, for instance, in the Lax-Wendroff theorem. We find 

Ij(6) - 7 3 (0) = (i , {vd t 6 + f(v)d x 9)dxdt. (3.25) 

To deal with 7^ (9), we recall that the flux <7jfc( n ) +e h( n )_i/2 depend on two arguments 
so satisfies 

f J£(n)+ e $(n)-l/2 = /( u Jfc(n)+ e *(n)) + 0(|Ujh (n) _ 1+e fc (n) - Mjh(„) +e fc(„)|) ^ 
= /(™q(*n)) + 0(|Wj£ (n )_ 1+e J( n ) - "jfc(„)+ e h(„)l)- 

Indeed, by construction, the point (n) + e^(n) — 1/2 is an end point of a cell 
achieving an extremum value: One has e^(n) G {0,1}, so J q (n) + e q (n) — 1/2 e 
{J*(ra) - 1/2, Jj(n) + 1/2}. Using (3.26), estimate (3.17), the Cauchy- Schwartz 
inequality, and finally Lebesgue convergence theorem, it is not hard to prove that 

I 2 h (0) -> 7 2 (0) = / f(w q (t))9(^ q (t), t)dt. (3.27) 

It remains to prove that 

1^(9)^1,(9) = - j R ^(t)w q (t)9(Mt),t)dt. (3.28) 
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We return to the definition of the modified path and define e q (n + 1) by 

J h q {n + 1 ) + e h q {n + 1 ) = jj(n) + ej(ra) + e h q {n + 1). 
Using only (3.21), one can prove 

E e J( n ) u "f(n) +e 5(„)^(„ )+e j (n) fc - - / *(f)0(V> ff (f),f)dt. (3.29) 

We claim that 

ej(n+l)e{ -1,0,1}. (3.30) 
Namely, using (2.11) and the definition of e q (n + 1), we have either 

ej(n) = 1 Jj(n + 1) = J q (n) or + 1) + 1, and e h An + 1) = 1 or 1 



or 

ej(ra) = 1 Jg (n + 1) = J^n) or j£(n + 1) + 1, and e h q (n + 1) = -1 or 0. 
The term l£(6) then can be rewritten in the form 

i*(o)= e u] +i ej +i h- Yl u l e 7 h 

j>J*(n)+e$(n) j> J* (n)+ e J (n) 

E E w 

j>Jh(n-l)+e£(n-l) i>J^(n)+eJ(n) 

so that 

= " E U ^(n)+6^(n)^J^(n)+eJ(n)^+ E U J, fc (n)+e5(n)-l^(n)+e5(n)-l' 1 - 
e*(n) = -l ej(n)=l 

Observe that, in the second sum above, J q (n) + e q (n) — 1 = (n — 1) + e q (n — 1), 
and consider the decomposition I^{9) = lii(0) + /^(fi) with 

^1,1 W = E ( U J£(n)+e5(n)-10Jf(n)+e£(n)-l ~~ U J£ (n)+ e J (n)# j£ (n)+ e J (n)) /l 
ej(n) = l 

and 

IlA e )= E U J^(n)+e5(n)^Jf(n)+e5(n)' 1 - E U J* (n)+e£ (n) ^ (n)+ e J (n) /l - 
e£(n)=l e£(n) = -l 
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On one hand, we have 

< Yl l^ h (n)+ e J(n)-l^(n)+^(n)-l ~ U Jh( n )+ e h( n )0jh (n)+e h (n) \h 
ej(n) = l 

<0(1) ^ |«Jf(n)+ e J(n)-l - u Jf(»)+ e J(»)l /l 
e*(»)=l 

+ 0(1) Yl \ B J£(n)+ e $(n)-l ~ d J£(n)+e%(n)\ h ' 
ej(n) = l 

and, in view of (3.17) and the smoothness property of 6, 

I} A {6)< 0(l)^ 1/2 (El«"f(n)-l -«"j(n)| 2 ) 1/2 + 0(l)^ 

n 

<0(1)(^ + ^/ 2 ) < 0(l)/i 1/2 , 

which implies 

7^(0) ^ as /i -> 0. (3.31) 
The expression for 7^2 (#) can be simplified, namely 

^1,2(0) = -E e !(") M "j(«)+ £ ;(»)^(™)+ £ ;(«)' i ' 

n 

Using (3.29), we find that 

7 i,2(^) h{0) as /i -> 0. (3.32) 

In view of (3.25), (3.27), and (3.28), we conclude that 

h(6) + 1 2 (6) + h(6) = 0. (3.33) 

Finally, using in (3.33) a sequence of test-functions 6>, whose supports shrink and 
concentrate on the curve ip qi the desired Rankine-Hugoniot relation (3.20) with the 
"+" sign follows at the limit. This completes the proof of Lemma 3.8. □ 

Proof of Lemma 3.9. For definiteness, we assume that u n Th/ >, is a maximum 
value and that: 

min(-Ujh (n) - Ujfc( n ) ±1 ) = "jfc(n) ~~ u "fc(n)-i- 

The other cases are treated similarly. To simplify the notation, set = Jq{n). By 
the uniform decay property (2.12), we have 

>«(«"*" «".-i) 2 - ( 3 - 34 ) 
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By the strong maximum principle (2.11), we have 



ut-i-ul-i<\K (3-35) 



From (3.35), we deduce 



«r-i<<-i+^«-«]f.-i) 



<-o(«Mi) 



< «£ - - min ((«£ - «£_!), («£ - ^.i) 2 ) 



n+l 



(3.36) 



Note that — _ 1 might be either < 1 or > 1 . Then we get 

^ - ^ min (« - K* - ^-i) 2 )- (3-37) 

On the other hand, (3.34) can be written in the form 

,f ^J.-a^-^f. (3.38) 

Moreover, using (2.11) again, we obtain: 

Hence, in view of (3.36): 

< u l - Z 3 ™ (K - K - ^-i) 2 )- (3-39) 

It follows from (3.37)-(3.39) that 

max(u]+_\ lU ]+\u]+ + \) <ul -/?min(« -<_!),« -<-i) 2 )> ( 3 -40) 
and thus 



?(n + l) < < " (K - « - «5l-l) 2 ). (3-41) 



since ^jh"( n+1 -) by definition achieves the maximum in (3.20). Finally, we have 
proved 

U W - «3f (Ui) ^ ^ min (« - K* - <-i) 2 )- (3-42) 

By assumption, it!? — u^-i < — . Thus (3.42) gives (3.17) after summation 
w.r.t. n. □ 
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Lemma 3.10. We have 

v(x—,t) > v(x+,t) in the set Qs(v). (3.43) 

□ 

The proof is based on the fact that one of the entropy criteria is satisfied on 
each side of a path: the Lax shock inequality in the non-increasing side, and the 
cell entropy inequality in the non-decreasing one. 

Proof. We claim that, along any path of extremum ip q , 

v(il> q (t)-,t) > v(if> q (t)+,t) for a.e. t E {0,T q ). (3.44) 

We use the notation introduced in the proof of Lemma 3.8. A new difficulty arises: 
several paths may accumulate in a region. Lemma 3.8 was concerned with the 
discrete conservation laws (2.1) which holds in both the non-increasing and non- 
decreasing regions. For the entropy consistency, we do not use the same criterion, 
and this complicates the proof. 

To begin with, consider a path ifj q and a point (ip q (to),to), that is supposed to 
be an "isolated" point of change of monotonicity, in the sense that: ip q -i(to) < 
iiq(to) < V , «+i(*o)- By continuity, these inequalities then hold with to replaced by 
any t lying in a small neighborhood of to- F° r definiteness, we also suppose that 
ip q is a path of minimum values. We later analyze the case that two or more paths 
of extrema accumulate in a neighborhood of (ip q (to), to). 

Let 9 be a non-negative test-function of the two variables (x, t) having its support 
included in a small neighborhood of (ip q (to),to)- We can always assume that u h 
is non-increasing on the left side of the curve ip q , and non-decreasing on the right 
side. Using the notation introduced in the proof of Lemma 3.8, we aim at passing 
to the limit in 

I h (0) = - E (U? +1 -U? + \(G] +1/2 -G]_ 1/2 ))9?h>0. 

Note that I h (9) is non-positive according to the cell entropy inequality (2.9) and 
since the u h, s are non-decreasing on the right side. 
Integrating by parts in I h {9) gives 

I h (Q) = ~ E (U™ +1 9] +1 - U?9™)h + G n jh {n)+e h {n) _ 1/2 9 n jH {n)+e h {n) _ 1/2 r 
+ £ U- +1 (9J +1 -9-)h + G] +1/2 (9- +l -9-)r 
= I?(9) + I 2 h (9)+I%(9). 

(3.45) 
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The passage to the limit in the term 1^(9) is a classical matter. The treatment of 
I 2 (9) and Ii(6) is similar to what was done to prove (3.27) and (3.28), respectively. 
Therefore we have 

!?(*)-> Ii(0) = - ^-(t)U(w q (t))9(^ q (t),t)dt. (3.46) 

/£(0)-/ 2 (0)= / F(w q (t))9(^ q (t),t)dt, (3.47) 
7r+ 

-> 7 3 (0) = // r , {U(v)d t 9 + F(v)d x 9)dxdt. (3.48) 

J J {x>ip q (t)\ 

It follows from (3.46)-(3.48) that 

h(9) + h(9) + I 3 (6) > 0. (3.49) 

Finally, using a sequence of test-functions whose supports shrink and concentrate 
on the curve ip q , we deduce from (3.49) that the entropy dissipation is non-positive 
along the path, i.e. 

-^(t)(U(v(i> q (t)+,t)) -U(w q (t))) + F(v^ q (t)+,t))-F(w q (t)) < 0. 

Combined with (3.20), this inequality is equivalent to 

v(%(t)+,t))<w q (t), 

which yields the desired inequality (3.44). 

Consider next the case that several paths accumulate in the neighborhood of 
(ip q (to),to). In view of (1.4), a finite number of paths only can accumulate at a 
given point. For definiteness we suppose that 

- the point (ip q (to), to) is a point of minimum values for u ; 

- the curves V'g+i an d coincide with ifj q in a neighborhood of to; 

- and we have iftq-i < ^ q and ip q +z < m a neighborhood of to. 
Suppose that, for instance, w q (t) < w q+ 2(t). The other cases are treated similarly. 
In this situation, we are going to prove that 

v(4> q (t)-,t) > w q {t) = w q+1 (t) > w q+2 (t) = v(i!> q (t)+,t), (3.50) 

at those points t near t where ip q (t) = ip q+ i(t) = ip q+2 (t). Of course, (3.41) is a 
much stronger statement than (3.34). By definition of the paths of extrema, we 
have 

w q {t) < w q+1 (t) and w q+1 (t) > w q+2 {t), 
v(ip q (t)-,t) > w q {t) and w q+2 (t) < v(ip q (t)+,t). 
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Lemma 3.8 shows that 

w q (t),w q+1 (t), andw q+1 (t) e {v{ip q {t)-,t),v{ip q {t)+,t)}. 
Thus, in order to get (3.50), it is sufficient to check the following two inequalities 

w q+2 (t)>v(if> q (t)+,t), (3.51) 

w q (t) > w q+1 (t). (3.52) 

On one hand, the argument used in the first part of the present proof applies 
directly to the path ip q +2(t) and the region located to the right of this curve, 
since ip q +2(t) < ip q +3(t) in a neighborhood of to- As a consequence, we obtain 
w q+2 (t) > v(ip q +2(t)+,t), which is exactly (3.51), since ifi q +2 = Vv 

On the other hand, to prove (3.52), let V q q+1 be the (small) region limited by 
the curves ip^ and ip q+1 for t belonging to a small neighborhood of to. Specifically, 
V q q+l is a set of indices of the form (j, n) defined along the lines of the proof of 
Lemma 3.8. In particular, both paths are modified according to estimate (3.17), as 
explained before. Consider 

I\0) = - J2 {U? +1 -U? + \(G] +1/2 -G]_ 1/2 ))e?h>0. (3.53) 
Note that I (9) is non-negative and 

i h {e) = - {u™ +1 ej +1 - u?e?)h 

I \ ^ an \ ^ rm an 

+ ^Jf (n)+eJ(»)-l/2 U J q -(n)+eh(n)-l/2 2^ U J* +1 (n)+e£ +1 (n)-l/2 ° (n)+e£ +1 (n)- 

n n 

+ E - + G i+i/2 w+i - e i> 

U,n)€V* q+1 

= ii{e) + % q {e)-iS, q+1 (e) + i$(e). 

Using the technique developed for the proof of Lemma 3.8, we get 
I}(0) - h{9) = -J r ^(t)U(w q (t))0(4&),t)dt 

+ J r ^^(t)U(w q+1 (t))e(i; q+1 (t),t)dt, 
lU e ) - W*) = / F{w q (t))6U; q (t),t)dt, 

i2, q+ M -> i 2>fl +i(0) = - / FMt^^iW.t)*, 

J 3 ft (0)^O. 
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It follows that 




(3.54) 



□ 



It is a classical matter to check that the initial condition (1.2) is satisfied by 
using (2.2) and the uniform BV bound. Since the function v(t) has bounded total 
variation, it admits left and right traces at each point and (1.5) and (1.6) are known 
to be equivalent at a point of discontinuity. Therefore the following result follows 
from Volpert's proof in [37]. 

Lemma 3.11. Supose v is a function of bounded variation and a weak solution to 
the conservation law (1.1) and satisfies the initial condition (1.2), and the inequal- 
ities (3.14), (3.15), and (3.34) where Q^v) U Q 2 {v) U Q 3 (v) = R x R + . Then v is 
the unique entropy solution to (1.1) -(1.2) . □ 

The proof of Theorem 2.3 is now complete. 

The limiting paths ip q associated with the scheme determine a decomposition 
of the plane into non-increasing/non-decreasing regions for the exact solution u. 
Such a decomposition is not unique, in general. Consider the decomposition found 
in the Appendix for the function v and the corresponding paths ip p . When v is 
not constant in any neighborhood of an extremum path <p p , the path is unique and 
must coincide with one of the path ip q . When v is constant in the neighborhood of 
a path <f p , then the path may be arbitrarily modified and it may happen that no 
limiting path ip q coincide with cp p . 



The purpose of this section is to apply Theorem 2.3 to van Leer's MUSCL scheme 
(for Monotone Upstream Scheme for Conservation Laws); cf. [23, 24]. This section 
also provides a proof of estimate (2.12) stated in Proposition 2.1, a new property 
of the Godunov scheme which does also hold for the MUSCL scheme. 

It is convenient to formulate (2.1) in terms of the incremental coefficients Cj^^ 2 
defined by 



4. Application to the MUSCL Scheme 




97+1/2 ~ /K) c -, n /("?)- ff-1/2 



(4.1) 





so that 



= u] + C+£ /2 (u? +1 - u?) + <7_l?/2(«"-i - «?)• 



(4.2) 



The numerical viscosity coefficient (Cf. Tadmor [35]) being defined by 



7 + 1/2 — °j + l/2 ' °j+l/2' 



(4.3) 



the viscous form of the scheme is 



< +1 = u n ,-\ (/K +1 ) -f(uU)) + ^W«£k -«?) - \QU/2 ( 



u 



n 
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Proposition 4.1. The scheme (2.1) satisfies the local maximum principle (2.11) 
provided 

q^ /2 >0, Cy- /2 >0, and Cfr» /2 + C7» /2 <± (4.4) 

□ . 

A sufficient condition for (4.4) to hold is 

Cf£ /2 >0 and QJ +1/2 <l/4. (4.5) 

Namely, if (4.5) holds, then < C+^ /2 < 1/4 and < C~£ /2 < 1/4, so that (4.4) 
is satisfied. In particular, the Godunov and Engquist-Osher schemes satisfy (4.4) 
under the CFL condition (2.4). When the numerical flux in independent of A, the 
second inequality in (4.5) is always satisfied provided A is small enough. 

The Lax-Friedrichs type scheme have a constant numerical viscosity Q™ +1 / 2 = Q- 
For the original Lax-Friedrichs scheme Q = 1. Proposition 4.1 applies provided 
Q < 1/4. Observe that the monotonicity property does fail when Q e (2/3,1]: 
take for instance / = and u™ = for all j ^ but Uq > 0. This initial data has 
one maximum point, and at the next time step 

U] +1 = + (1"QK +f<-l- 

admits two maximum points j = — 1 and j ; = 1. A related observation was made 
by Tadmor in [35]: for Q < 1/2, better properties can be obtained for the Lax- 
Friedrichs scheme. 

Proof of Proposition 4.1. Inequalities (2.11) can be written in terms of the 
incremental coefficients, namely 

\ min *» 5U) < ^i/ 2 ^i + " ^i/ 2 " ^-1/2)^ + ^-1/2^1 

<imax(^ +1 ,^,^_i) 

(4.6) 

with 

If (4.4) holds, then 

2 c;+? /2 *? + i + (1 - 2 ^1/2 - 2 ^7-1/2)^ + 2 c ;-i/2 5 u) 

is a convex combination of the Sj's. So (4.6) and therefore (2.11) follows. □ 

We now introduce the van Leer's scheme, composed of a reconstruction step 
based on the min-mod limitor and a resolution step based on the Godunov solver. 
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We use the notation introduced in Section 2. For simplicity in the presentation, we 
normalize the flux to satisfy /(0) = f'(0) = 0. From the approximation [u 7 - } at 
the time t = t n , we construct a piecewise afRne function 

Uj(x) = u 7 - + sl](x - Xj)/h for x e (x j _ 1 / 2 , x j+1 / 2 ), (4.7) 

where the slope s 7 - is 

s" = minmod(wj — itj_i, (%+i — %-i)/2, %+i — %) (4-8) 

with 

min(a, 6, c) if a > 0, 6 > 0, and c > 0, 
minmod(a, 6, c) = ^ max(a, 6, c) if a < 0, b < 0, and c < 0. (4-9) 

in all other cases. 

We introduce the notation 

mJ +1/2 _ = + ^ n /2 and u] +1/2+ = u] +1 - s] +1 /2. 

Then the solution is up-dated with (2.1) where the numerical flux is defined de- 
pending upon the values of the reconstruction at the interfaces. 

(1) If either < u™ < or < Uj +1 > u™, then the numerical flux is 
defined by using the characteristic line traced backward from the point 
(xj+i/2,i n +i/2)- Since the latter has a positive slope, we set 

9?+i/2 = /K+i/2-) + Z'K+i/2-) (v - «J +1/2 _) (4.10a) 
with v solving 

u] +1/2 _ = v + ±f'(v)s]. (4.10b) 

(2) If either < Uj +1 < 0, or u™ +1 / < ulf < 0, then the backward characteris- 
tic has a negative slope and we set 

9?+i/2 = /«+i/2 + ) + f'W +1/2+ ) (v - u] +1/2+ ) (4.11a) 
with v solving 

«?+i/2 + = « + ^/ / W+i- (4-llb) 

(3) In all other cases we set 

9j+i/2 = /(0) = 0. (4.12) 

Equations (4.10b) and (4.11b) can be solved explicitly for the Burgers equation 
since then f'(u) = u is linear. Observe that the scheme reduces to first order at 
sonic points and extrema. 

The main result of this section is: 
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Theorem 4.2. For A small enough, the MUSCL method defined by (4.7) -(4. 12) 
is a generalized monotone scheme in the sense of Definition 2.2. When (1.3)-(1.4) 
hold, the scheme converges in the strong L 1 topology to the unique entropy solution 
o/(l.l)-(1.2). ' □ 

It would be interesting to extend Theorem 4.2 to higher-order methods such as 
the Woodward-Collela's P.P.M. scheme. 

Proof of Theorem 4.2. We have to check that the scheme satisfies the three 
conditions in Definition 2.2. We always assume that A is, at least, less or equal to 
1/4. 

Step 1: Local maximum principle. 

Estimate (2.11) is easily obtained by applying Proposition 4.1 and relying on the 
convexity of the flux function /. We omit the details. 

Step 2: Cell entropy inequality. 

Consider a region where the sequence {w™} is non-decreasing. We will use the 
entropy pair (U, F) with U(u) = u 2 /2 and F'(u) = uf'(u). Define the numerical 
entropy flux by 

G?+l/2 = FW+1/2-) + ^K+l/2-) (/(«) - /K+i/2-))> (4.13i) 
G?+i/2 = *K+i/ 2+ ) + U'(u] +1/2+ ) (f(v) - /K +1/2+ )), (4.13ii) 

and 

G7+1/2 = ( 4 - 13iii ) 

in Cases (1), (2), and (3), respectively. Inequality (2.9) is checked by direct calcu- 
lation, for A small enough. Observe that Case (3) is obvious since our scheme then 
reduces to a first order, entropy consistent scheme. 

For definiteness we treat Case (1), i.e. /' > for the values of u under consid- 
eration. We view the left hand side of (2.9) as a function of w = Uj_i/ 2 -, u = U™, 
v = Uj + i/ 2 - and the value w defined as 

w = w + ^f'(w)t, 

where t stands for the slope in the cell j — 1. Introduce also v by 

v = v + ^f'(v)s 

with s = 2(v — u). Since the approximate solution is non-increasing, we have 
w < w < 2u — v < u < v < v. Set 

Q(w, to, u, v; A) 

= U{u) - U{u) +X[F(v) + U'(v)(f(v) - /(«)) - F(w) - U'(w)(f(w) - f(w))], 
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and 

u = u- X[f(v) + f'(v)(v -v)- f(w) - f'(w)(w - w)] . 
Observe that 

d^Q = U'(u)Xf'(w) - XU'(w) f'(w), 

and 

dlfl = U"(u)X 2 f'(w) - XU'(w) f"(w) < -CX \w\ 
for w > and A small enough. Therefore O is a concave function in w and 

Q(u), -u, v; A) < fi(iy, w, u, v; A) — (w — w, -u, v; A) — CA|w| \w — w\ 2 . 

But 

w, u, v; A) = Xf'(w)(U'(u) — U'(w)) 

= U"(0 A/» [u-w- X(f(v) + f'(v)(v -v)- f(w))} 

< CX \u — w\ 

for some £ > and A small enough. This proves that 

Q(u), w 7 u, v; A) < 0(w, w, u, v; A) — CX(\w — w\\u — w\ + \w\ \w — w| 2 ) 

< Q(w, w, u, v; A), 

and we now simply use the notation -u, v; A). 

Taylor expanding O with respect to A shows that the dominant term is the first 
order coefficient in A given by 

Ox (w, u, v) = -U'(u) (f(v) - f(w)) + F(v) - F(w), 

in which w <2u — v < u < v. Since 

w fiiKu,v) = (U'(u) - U'{w)) f'(w) > C\w\ \u-w\, 

we have 

flu — v 

tii(w 7 u,v) < Oi(2w — v,u,v) — C \w\ I (u — z)dz, 

J w 

SO 

£l\(w, u, v) < &\{2u — v,u,v) — C \ w\ \2u — v — w\ \u — w\. 
It remains to study Cli(2u — v, u, v) = Qi(u, v) with u < v. We find 

d v Qi(u,v) = (U'(v)-U'(u))f'(v)+ (U'(2u-v) - U'{u))f'{2u- v) 
= (v-u) (/» - f\2u - v)) <-C\u- v\ 2 . 

It follows that Q(u, v) is a non-increasing function of v for all v > u, and since it 
vanishes for v = u, 

Cli(u,v) < -C \u-v\ 3 . 
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This proves that the first order term in A in the expansion of the function O is 
negative. 

The same arguments are now applied to the function 0(A) directly. We have 

d w tt ( w, u, v; A) = (U'(u) - U'(w)) Xf'(w) 

= \f(w)U"(0(u -w- X(f(v) + f(v)(v -v)- f(w))) 
= Xf(w)U"(£)(u -w- XO(l)(u - w)) 
> AC \w\ \u — w\ . 

Therefore, 

Q(w, u, v; A) < Vt{2u — v, u, v; A) — AC' \w\ \2u — v — w\\u — w\. 
Denote Q(u, v; A) = Q(2-u — v, u, v; A) 

Cl(u, v; A) = U(u) - U(u) + A [F(v) + U'(v)(f(v) - f(v)) - F(2u - v)] , 

where 

u = u- X[f(v) + f'(v)(v -v)- f(2u - v)]. 
We easily compute that 

^ fl( ^ ;A) = f'(v) (U'(v) - U'(u)) + f'(2u - v) (U'(2u -v)- U'{uj) + XA(u, v; A) 
A 

with 

\A(u,v; A) | < C\u- v\ 2 . 

This establishes the desired conclusion for A small enough. 

Step 3: Quadratic decay property. 

Near a local extremum, the MUSCL scheme essentially reduces to the Godunov 
scheme. So it is enough to check the quadratic decay property (2.12) for the Go- 
dunov scheme. This can be done from the explicit formula (2.6). 

The simplest situation is obtained with the Godunov scheme and when /' has a 
sign, say is positive. Assume u™ is a local maximum. We have 



1 = 



A(/(«?) -/(«?_!)), 

thus 



a(/K)-/K_x)) 

> A/'K_J («» - u? J + A(inf /72) («» - u") 2 



>A(inf/"/2) min(^-^ ±1 ) 2 
This establishes (2.11) when f > 
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It remains to treat the sonic case where /' has no definite sign. We will rely on 
the following technical remark. Given two points such that 

u- < < u+, /(«_) = /(«+), 

there exist c\,ci > (independent of u±) such that 

Cl \U-\ < \u + \ < C2 \U-\. 

Consider the case u 1 J_ 1 < < u™, and use Osher's formula for the Riemann 
problem, we have 



it? - -u? +1 = A ( max / - min /) 
= A (, max f-/(0)) 
= A(inf/72)( | max /) 2 



> c|u? +1 -u?| 2 . 



Consider next the case < < tf, then 

u?-u] +1 = \( max f-f(u])) 

> A(inf /' 72) min(|^ +1 - ^ n | 2 , - ^.J 2 ) 

This completes the proof of Theorem 4.2. □ 
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Appendix: Monotonicity Property. 

In this appendix we brieffly discuss the monotonicity property together with more 
basic properties of entropy solutions to conservation laws, which go back to Kruzkov 
[19] and Volpert [37]. In the paper by Keyfitz, [18], somewhat simpler proofs are 
available for piecewise Lipschitz continuous solutions. We are interested in the 
local versions of the properties, i.e. formulated in domains limited by characteristic 
curves. To cope with discontinuous solutions, we use the concept of generalized 
characteristic curves introduced for ordinary differential equations by Filippov [9] 
and developed in the context of conservation laws by Dafermos; see e.g. [7] and 
the references therein. We recall that, through any point (xo,to), there exists a 
funnel of forward and backward generalized characteristic curves, which fill up a 
domain (£ m (t) < x < £ M (t)}. Here £ m (respectively £ M ) is called the minimal 
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(resp. maximal) characteristic curve originating at (xo,to). It is known [9, 7] that 
a characteristic, say £, is Lipschitz continuous and for almost every time t > 

f rwm,t)) a u-(t) = u+ (t) = u(at),t), 

* (<) i^^ (Ai) 

where u± = it(£(£)±,t). For our purposes, / is strictly convex and there is a 
unique forward characteristic issued from (xq), to and there is no need to distinguish 
between the minimal characteristic and the maximal one, with the exception of 
those points where to = and uo has an increasing jump at xq] cf. [7]. 

Solutions u to (1.1) are Lipschitz continuous in time with values in L 1 and, for 
all times t, u(t) has bounded total variation in x. 

The following properties follows from [37, 19] and the technique of generalized 
characteristic in [9, 7]. 

Proposition A.l. Let u be the entropy solution to (1.1) -(1.2) . Given x\ and X2 
with x\ < X2, consider the maximal forward characteristic £" (t) issued from (0, x\) 
and the minimal forward characteristic ££(£) from (0, £2)- For all times t > 0, u 
satisfies 

(1) the local maximum principle for all t < s and y G (£1 (s), £2 ( s )) ; 

c «^ inf eunA^'Q} ~ u ( s ' y ^ ~ sup ( A - 2 ) 

(2) the local L 1 contraction property; 

d f min($(t),&(t)) 

— \u(x,t) -v(x,t)\dx <0, (A.3) 

(3) the local order preserving property; 

if u (x) < v (x) for all x G (max(£?(0), £(0)), min(£ 2 w (0), ££(0))), 

then u(x,t) < v(x,t) for all x G (max(fj*(t), £!"(£)) < x < min (££(£), ££(£))), 

' (A.4) 

(4) and the local TVD property; 

^TVg$«(f)<0. (A.5) 

where in (1) and (3) the function v denotes the solution to (1.1) -(1.2) with uq 
replaced by a function vq G BV(R), and the curves £1 (t) and ££(£) are defined in 
the obvious way. □ 

In (A.5), TY b a (w) denotes the total variation of a function w : (a, b) — > R. The 
derivatives in (A.3) and (A.5) are to be understood in the distributional sense. At 
least with x\ = £" = —00 and X2 = £2 = °°5 (A-5) is a direct consequence of (A.3) 
and the invariance by translation of the solution-operator for (1.1). 
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The following proposition concerns the monotonicity property, which is a refine- 
ment to the statement that the solution-operator is monotonicity preserving , i.e. 
satisfies 

if uo is monotone, then u(t) is monotone for all t > 0. (A. 6) 

Namely (A. 6) follows from (A. 3) by taking x\ = = — oo and x 2 = = oo and 
vo(x) = uo(x + y) for positive or negative values of y. 
It is convenient now to assume that (1.4) is satisfied. 

Proposition A. 2. (Monotonicity property) Suppose that the initial condition u$ 
has a locally finite number of local extrema. There exist (Lipschitz continuous) 
generalized characteristic curves tp q : [0, oo) — > R -the index q describing a subset 
E(uq) of consecutive integers- such that 

¥q < <Pq+l, (A.7) 

there is only a finite number of such curves in each compact set (A. 8) 
and, for all t > and all relevant values of p, 

u(t) is non- decreasing for x G (<f2p(t), <f2p+i(t)), 

u(t) is non-increasing for x & (ip2 P -i(t), (p2 P (t)), ^ 

and, as long as <fi2p(t) ^ (p2 P +i(t), 

u(if2p(t)+,t) is non- decreasing, 

v . . . (A.10) 

u(ip2p+i(t)—,t) is non-increasing. 

□ 

The paths </?2 P and v^p+i are called path of local maximum and path of local 
minimum for u, respectively. By convention, ip q = — oo for q > ma,xE(uo) and 
ip q = +oo for q < mmE(uo). 

If two initially distinct paths cross at a later time, then from that time they will 
coincide thanks to the uniqueness property for forward characteristics. Note also 
that such paths need no longer be paths of local extremum in a strict sense, but 
arbitrary characteristics , even though (A. 9) would still hold. On the other hand, 
when simultaneously uq has a decreasing jump at a point xq, that uq(xq— ) is a local 
maximum, and u(xq+) is a local minimum, then two equal paths originate from 
(0,xo), one being a path of minimum and the other a path of maximum. When 
simultaneously uq has an increasing jump at xo, that uo(xq— ) is a local minimum, 
and that uo(xo+) a local maximum, then two distinct paths originate from (0, xo). 

Proposition A. 2 is a classical matter in the literature although no specific refer- 
ence seems available. Cf. however Harten [13] and Tadmor [35] where the ideas are 
developed. 
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Proof of Proposition A. 2. First of all, the points <p q (0) and the set E(u ) are 
defined from the initial condition uq in an obvious way so that the conditions (A. 7)- 
(A.10) hold true at time t = 0. Let us define f2p(t) to be the maximal forward 
characteristic issued from <p2 P (0). Similarly, let ip2 P +i(t) be the minimal forward 
characteristics issued from </?2p+i(0). Indeed one need to distinguish between mini- 
mal and maximal characteristics only in the case of an initially increasing jump. It 
may happen that both a path of minimum and a path of maximum may originate 
from such a point of increasing jump. 

Property (A. 7) is an immediate consequence of the uniqueness property of the 
forward characteristic. Condition (A. 8) follows from the property of propagation 
with finite speed satisfied by solutions to (1.1) and the fact that the initial data has 
a locally finite number of local extremum. Indeed (A.l) yields a uniform bound for 
the slopes of the characteristics. 

In order to establish (A. 9) and (A. 10), we first suppose that uq does not admit 
increasing jumps. Consider an interval of the form (<f2 P (t), <fi2 P +i(t)) for those 
values of t when this interval is not empty. Note first that, taking £J* = ip2 P and 
£2 = f^p+i, the local maximum principle (A. 2) implies in particular that 

u(ip 2 p(t) + ,t) > U (lfi2p(0)+), 
U(<fi2p+l(t)-,t) < U (<P2p+l(0)-). 

Let w be the entropy solution to (1.1) associated with the initial condition 

«o(</?2p(0)+) if X < <fi2p(0), 

w(x, 0) = w (x) = u (x) if v?2 P (0) < x < <f2p+i(0), (A.12) 

UO (^2p+l (0) -) if X > if2 P +l (0) . 

The data too is non-decreasing and, in view of (A. 6), 

the solution w is non-decreasing for all times. (A. 13) 

Let •02 P and V^p+i be the forward characteristics associated with w and issued 
<f2 P (0) and </?2 P +i(0) at time t = 0, respectively. Observe that the maximum forward 
and the minimum forward curves coincide since by construction wq is continuous 
at </?2p(0) and <p2 P +i(fy- Note in passing that the function w satisfies: 

I U {<P2 P +l{0)-) llX> 1p2p+l{t). 

Using (A. 11) and (A.14) and the fact that /'(.) is increasing, one gets 

< fWv2 P+ l(t)-,t)) < f(u (v2 P+ l(0)-)) < d ^(t), 



which implies 



V2p+l(t) < V>2p+l(£)- 
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Similarly 



Using the L 1 



contraction principle (A. 3), it follows that 



u = w 



for <p2 P (t) < x < (p 2 p+i(t), 



and, in view of (A. 13), the function u is non-decreasing and (A. 9) holds. Using 
(A. 9) and the local maximum principle (A. 2) finally provides (A. 10). The proof 
is complete in the case of an interval of the form (<p2 P (t), <f2 P +i(t)). An interval 
(<f2 P -i(t), <f2 P (t)) can be treated in a similar fashion. 

It remains to consider increasing jumps in uq. That situation can be treated by 
using the following property. Suppose uq has an increasing jump at a point xq and 
let (p m (t) and (p M (t) be the minimal and maximal forward curves from xq. It is 
known that at least for small times the function u(t) coincides with the rarefaction 
wave connecting the values u (x ±) in the interval (y? m (t), y? M (t)). 

This completes the proof of Proposition A. 2. □ 



